Magnetization and dimerization profiles of the cut two-leg spin ladder 
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^ , The physical properties of the edge states of the cut two-leg spin ladder are investigated by means 

Q ' of the bosonization approach. By carefully treating boundary conditions, we derive the existence of 

' spin-1/2 edge states in the spin ladder with a ferromagnetic rung exchange and for the open spin-1 

Heisenberg chain. In contrast, such states are absent in the antiferromagnetic rung coupling case. 
CO ' The approach, based on a mapping onto decoupled semi-infinite off-critical Ising models, allows us 

(N ■ to compute several physical quantities of interest. In particular, we determine the magnetization 

and dimerization profiles of the cut two-leg spin ladder and of the open biquadratic spin-1 chain in 
the vicinity of the SU(2)2 WZNW critical point. 

D 

i", PACS numbers: 75.10.Jm 75.30.Hx 75.50.Ee 71.10.Pm 

. I. INTRODUCTION 



The infuence of impurities and imperfections on the behavior of low-dimcnsional strongly correlated systems has 
attracted considerable attention in recent years. The introduction of static non-magnetic impurities, like Zn or Li, 
at the location of the magnetic ions is a sensitive probe of the correlations that develop in these magnetic systems. 
A particularly striking exemple is the observationrnof fractional spin-1/2 edge states in the Haldane gapti spin-1 
compound NENP cut by non-magnetic impurities.H These spin-1/2 degrees of freedom are associated with static 
staggered raoments close to the chain ends which are revealed unambiguously in the NMR profile of the Mg-doped 
Y2BaNi05.u This effect can be viewed as a local restoration of antiferromagnetism by impurities. In fact, a transition 
. to an antiferppmagnetic state induced by the local moments has been observed recently in the Haldane gap compound 
' PbNi2V208.Q Such a local enhancement of antiferromametism induced by non-magnetic impurities is a rather general 
T-H , phenomenon in gapful quasi-one dimensional systems.u The spin ladder material SrCu203 lightly doped with Zn 
' impurities exhibits a Curie-Uke behavior at low temperature which has been explained by the unpaired free spins in 
' the vicinity of the impurity.B Further, it has been shown by NMR measurements that a staggered magnetization is 
induced along the leg by very smaU-|tt).25%) concentrations of impurities. lJeI At low temperature, the induced moments 
get frozen, leading to NeeLprder.Elflla Similar effects have been observed in-the ladder compound Cu2(C5Hi2N2)2Cl4 
"j^ ■ doped with Zn impuritieslij and in the spin-Peierls gap material CuGe03.Ej 

■ A simple explanation of the existence of fijee spin-1/2 moments at the ends of a broken spin-1 chain can be obtained 
I from the Valence Bond Solid (VBS) modcO where each S=l spin is viewed as two S=l/2 spins in the symmetric 
''O ' triplet state. In this model, the singlet ground state of a chain with periodic boundary is described by two valence 
Ph , bonds originating from each site to form singlets with adjacent neighbors. If the chain is broken (i.e. open boundary 
O ■ conditions are considered), unpaired bonds are left at each end of the chain resulting in two free S=l/2 objects at the 
boundaries and a fourfold ground state degeneracy. This VBS picture proiddes a good and intuitive description of the 
ground state of the spin-1 chain. In particular, the exact diagonalizationlij of finite open samples with even number 
of sites has shown that the ground state is a singlet and the existence of an exponentially low-lying triplet state in 
the Haldane gap. This leads to a fourfold ground state degeneracy in the thermodynamic limit. Such degeneracy 
' can also be interpreted as the consequence of a spontaneously broken hidden Z2 x Z2 symmetrylia associated arith the 
formation of the Haldane gap. A string order parameter has been introduced to reveal this hidden symmetryEflllj More 
generally, it is expected within the VBS land non-linear sigma model approaches that the integer spin-S Heisenberg 
chain has spin-S/2 chain-end excitations.cZI 

The physical properties of the S=l/2 chain-bpundary excitations in the open spin-1 chain have been investigated in 
details in a^tuiaiijUira-M Carlo (QMC) studjil3 and by means of the density-matrix renormalization group (DMRG) 
approach. E2lE3'ElrE3c3 Recently, the magnetization profiles at finite temperatures and fields have been determined 
using continuous time QMC techniques to reconstruct the experimentally measured NMR spectrum of the Mg-doped 
Y2BaNi05.E3 The properties of the edge states of a more general model, the bilinear biquadratic spin-1 Heisenberg 
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chain, defined by the Hamiltonian: 



W = jV[S,-S,+i+/3(S, -S.+O^^], (1) 



have also been considered. For (3 — 1/3. (the so-called AKLT point) the VBS state turns out to be the exact ground 
state of the bilinear-biquadratic model.Ed For /3 = — 1, the Hamiltonian (|l|) has a critical point separating the Haldaae 
phase (—1 </?<!) from a dimerized phase {(3 < —1). At this (/? — —1) critical pppt, the model is integrablecS 
and belongs to the SU(2)2 Wess-Zumino-Novikov-Witten (WZNW) universality class.EZl The authors of Ref. have 
established that the edge states are present through the whole Haldane phase and disappear as soon as the (3 = —1 
critical point is reached. 

In this paper, we shall investigate the physical properties of the S=l/2 chain-end excitations ofrthe semi-infinite 
(or cut) two-leg spin ladder and of the open spin-1 chain by means of the bosonization method.c3 In the strong 
ferromagnetic rung limit, this two-leg ladder model is equivalent to the open spin-1 Heisenberg chain with the two 
spins on the rung forming an .effective S~l local moment. Since this strong coupling limit is smoothly connected to 
the weak coupling one,E3iElH the approach provides a simple way to extract the low-energy properties of the open 
spin-1 chain. In_Lhis respect, it gives an alternative derivation of-the existence of the spin-1/2 edge states predicted by 
the VBS theoryllj and the Schwinger-boson mean-field analysis.EZi Futhermore, we shall be able to calculate explicitly 
the physical propertieSj-of the open spin-1 chain such as the magnetization profile or the NMR relaxation rate. To 
this end, the mapping^ of the low-energy Hamiltonian of a weakly coupled two-leg ladder onto off-critical two- 
dimensional Ising models will he exploited to derive the chain-boundary excitations as it has been done in the study 
of disordered spin-1/2 ladders.L3 By paying a careful attention on the boundary conditions on the fields that occur in 
the continuum limit, the staggered magnetization profiles of the model can be determined using exact results of semi- 
infinite one-dimensional quantum Ising model. The results strongly depend on the sign of the interchain coupling and 
for an antiferromagnetic rung exchange, najmagnetic chain-end excitations are found. However, a weak dimerization, 
induced by the presence of the boundary,t3tirE3 exists for all signs of the interchain interaction and can be computed 
by this mapping onto semi-infinite Ising models. Finally, the infiuence of a strong external magnetic field fixing the 
spins at the edge can be investigated by a similar approach. 

The rest of the paper is organized as follows. In Section ||, the low-energy Hamiltonian of the cut two-leg spin 
ladder is mapped onto an 0(3) x Z2 symmetric theory of four massive Majorana fermions with suitable boundary 



conditions. The nature of the edge states that occur in the problem is then discussed in Section III where the uniform 



magnetization profile and the NMR relaxation rate are computed for a ferromagnetic interchain interaction. Section 



IV presents the calculation of the staggered magnetization and dimerization profiles of the model by exploiting the 
mapping onto semi-infinite off-critical quantum Ising models. The effect of a strong applied magnetic field fi xing the 
spins at the boundary is investigated in Section Finally, our concluding remarks are presented in Section VI and 
the paper is supplied with four Appendices which provide some technical details used in this work. 

II. DERIVATION OF THE EFFECTIVE LOW-ENERGY HAMILTONIAN 

In this section, we apply the bosonization approach to the semi-infinite two-leg spin ladder described by the 
Hamiltonian 



n=0 



P=l,2 



(2) 



where S„^p is a spin-1/2 operator at site n on chain p {p = 1,2) and we consider an antiferromagnetic inchain 
interaction Jy > 0. The bosonization method will be applied to the Hamiltonian (||) in the regime |Jl| ^ J|| and 
with a suitable redefinition of the effective coupling constants it captures the physical properties of the mOjdd fnr. 
arbitrary J± since there is a continuity between the weak and strong coupling limits in the two-leg spin ladder Ouj'EiJ 
In particular, local S=l spins are formed in each rung of the ladder in the strong ferromagnetic interchain coupling 
limit (Jj_ < 0) so that the approach provides, in turn, a way to investigate the physical properties of the broken spin-1 
Heisenberg chain. 

A. Bosonization of the open two-leg spin ladder 



Let us first consider the decoupling limit {J± = 0) where the system reduces to two independent spin-1/2 Heisenberg 
chains with open boundary conditions. The low-energy properties of the latter model can be still determined by 
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means of the bosonization method f H i As it is described in the Appendix A, starting from the underlying 

Hubbard model, the bosonized Hamiltonian for the open spin-1/2 Heisenberg chain with index p = 1,2 reads as 
follows neglecting the marginally irrelevant term: 



^7'' Jo 



(3) 



where v is the spin velocity and Hp is the momentum operator conjugate to the bosonic field <i>p. The boundary 



condition on the fields $p reads 



% (0) - 0, (4) 



which corresponds to a Dirichlct boundary condition. In the continuum limit, the effective spin density Sp(x) separates 
into an uniform and staggered parts: 

Sp [x) = JpR (x) + 3pL (x) + (-1)"=/" Hp (x) , (5) 

a being the lattice spacing. As shown in the Appendix A, the bosonized description for the uniform spin density is 
given by 



1 

:\/2*pR 



JpR,L - ^ dx'^pR^L 



27ra 

^pR,L being the chiral components of the bosonic field <i>p: $p = ($p7?, + <I>pl)/2. The staggered part of the spin 
density (|^) can be expressed in terms of $p and its dual field Qp-. 



A r 

TTfl 



IS (%/2 Gp) , - sin ( V2 Op) , - sin (^2 $p) ] , (7) 



where A is a constant stemming from the underlying charge degrees of freedom that have been integrated out. 

In the weak coupling regime |Jl| <C J\\, the continuum limit of the Hamiltonian (H) can then be derived using all 
these results. To this end, we introduce the symmetric and antisymmetric combinations of the bosonic fields: 

71 ' ^ ^ 

so that the leading part of thejBjftfpiltonian (||) that imposes the strong coupling behavior of the system decomposes 
into two commuting parts Tl± 



334 



n~n+ + n- [n+,n-]^o, (9) 

with the following bosonized expressions 

n+ = — dx [{TTll+f + {dx<P+f] - -V / dx COS 2$+ 
27r Jo 27r^a Jq 

H-^— da; [(7rn_)2 + (a^$_)2] + _L_ / ^fa; cos 2$_ + / da; cos 26- , (10) 

2n Jq 2TT^a Jq tt^o Jq 

where the boundary conditions on the bosonic fields are of Dirichlet type: 

$±(0) = 0. (11) 

In this derivation of the low-energy theory, one should note that we have only taken into account the most relevant 
perturbation that appears in the continuum limit of the spin ladder. In particular, we have discarded the marginal 
contribution that stems from the uniform pieces of the spin densities (|^). We shall later comment on its main effect 
when the Hamiltonian (|^) will be refermionized. 
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B. Refermionization 



The next step of the approach is to observe that the scahng dimension of the interacting part in Ti± is equal to 
one. The bosonic fields ar&^recisely at the free-fermion point where the cosine terms in Eq. (nG) can be expressed in 
terms of massive fermionsH3 To this end, we first introduce the left and right bosonic fields corresponding to 

$± = ^ (*±L + $±7?) 

e± - ^(*±L-$±fl,). (12) 



These chiral fields are no longer independent due to the existence of the boundary condition ( ]f f| ) and one has: 

$±^(0) = -$±l(0). (13) 
The refermionization of the Hamiltonian ( p^ ) can then be obtained through the bosonization formula 

V±R = ■ 



V 27ro 

= (14) 

where k± are Klein factors that obey the anticommutation relation — to ensure the anticommutation 

between the fermion fields with different channel index ±. The anticommutation between ijj±R and '0±l results from 
[^±B.ix),^±L{y)] = —ITT which stems from the Dirichlet boundary condition ( pT| ) as described in the Appendix A. 
The boundary conditions on the fermionic fields can be deduced from Eq. (p^ ) 

^±R (0) = V±L (0) . (15) 

The cosine terms of Eq. ( po|) can then be refermionized using the identification ( [T^ as well as the commutation 
relation [^±r{x),^±l{x)] = —iTT,x > 0: 

cos2$± = -iTTfl (^■0±i?V'±L - V'liV'ifl) 

cos2e± = JTi-a (^ViiiV^ii - V'±LV'±fl) • (16) 

The Hamiltonians 7i± of Eq. (^0|) can thus be expressed in terms of the fermion fields 

f°° I \ f°° I 

H+ = -iv J dx i^tpljidxip+R - i^l^d^^+Lj + ' J dx (^^^^-0+^ - -0^^?/^+^ 

f°° ( \ iJ±X^ f°° ( \ 
U- = -iv J dx i^ipii^dxi^^R- i^'l^dxij^Lj J dx \^^_riI^-l - '4>^-L'^^Rj 

rfo: (^l^Vli - V'-LV'-fl) , (17) 

with the boundary conditions ( [l5| ) for the fermion fields. 

At this point, it is convenient to introduce four Majorana (real) fermions from the Dirac ones (^4|) 

, ^R.L + ^^R,L 

i'-n.L = (18) 



This identification together with the correspondences (|6|,|14|) enable us to derive the refermionization of the uniform 
part of the spin densities (||) 



ra ja t ja ahc 
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^R,L = JlR,L ~ ^2fl,L = ^^R,L^R,L' (1^) 
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where we have fixed the product of the Klein factors that appear in Eq. (^_4|) to i to obtain Eq. (|19|). In 

fact, this identification ( pi)|) is nothing but the faitfifu l r ji pp resentation of two independent SU(2)i Kac-Moody currents 
•fpR,L,P = 1, 2 in terms of four Majorana fermions.uSEl 

With the above results, the Hamiltonian (^) can be rewritten with the four Majorana fermions and be separated 
into two commuting (triplet and singlet) pieces 

n = nt + ns, (20) 

with 

^ -^0 a=l -^0 

-Hs = r dx {end.a ~ iid^ii) - r dx aa. (21) 

^ Jo Jo 



The boundary conditions on the Majorana fermions are obtained from the constraint (y_5|) and the definition (18) 



eS(0)=?£(0), a = 0,...,3. (22) 

Moreover, the marginal interchain perturbation that we have so far neglected can be expressed in terms of the 
Majorana fermions using the correspondence (|l^). As shown in Ref. ^2[ the resulting contribution leads to a velocity 
anisotropy and a mass-renormalization in the singlet and triplet sectors so that the low-energy Hamiltonian ( |20| ) takes 
now the following form 

Hs = r dx {^nd^eR - ild.eL) - r dx ^r^l, (23) 

^ Jo Jo 

where uit > and rn^ < (respectively mt < and rris > 0) for a ferromagnetic (respectively antiferromagnetic) 
interchain coupling and in particular in the weak coupling case | Jl| <C J\\ one has the identification from Eq. (pi]): 
nit = — J_LA^/27r and nis — 3J±X'^/2n. 

We thus observe that in the low-energy limit the initial Hamiltonian (^) of the cut two-leg spin ladder is mapped 
onto a model of four free massive Majorana fermions with the boundary conditions (p^. In the strong ferromagnetic 
rung limit —J± S> Jy, the singlet excitation described by the Majorana fermion ^ are frozen (|m<;| — > oo) so that 
the low-energy properties of the model are governed by the triplet magnetic excitations corresponding to the fields 
^rlt^ = 1,2,3. In this strong ferromagnetic rung limit, we expect the system to be equivalent to a broken spin-1 
chain. Indeed, it was shown by TsvelikEJby perturbing around the SU(2)2 WZNW critical pointH of the biquadratic 
spin-1 chain, which is described by three massless Majorana fermions, that the low-energy properties of a gapped 
spin-1 chain could be obtained from a triplet of massive Majorana fermions. Furthermore, it can be seen easily that 
the boundary conditions (|2^ ) imply that the SU(2)2 currents obey /^(O) = /£(0) which means that there is no spin 
current flowing across the boundary. Therefore, an open biquadratic spin-1 chain is described by a triplet of massive 
Majorana fermions: 

/•oo 3 ^oo 3 

Ut = --^ dx Y.&-^R-^ld-CL)-^mt / dxY, (24) 

^ -^0 a=l -^0 a=l 

with the boundary conditions ( p2| ) and the Haldane (respectively dimerized) phase is characterized by a positive 
(respectively negative) triplet mass mt- 



III. S=l/2 CHAIN-BOUNDARY EXCITATIONS 

In this section, the nature of the edge states of the cut two-leg spin ladder and open spin-1 chain are investigated 
using the low-energy description ( p3| ) of the model in terms of four Majorana fermions with boundary conditions. In 
particular, physical quantities such as the uniform component of the magnetization profile and the NMR relaxation 
rate will be computed within this approach. The calculation of the staggered magnetization near the edge will be 
presented in the next section since it involves quantities that are nonlocal in terms of the Majorana fermions. 
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A. Localized Majorana fermion state 



The special structure of the low-energy Hamiltonian (^3|) together with the constraint (g2|) lead us to consider a 
single massive Majorana fermion Hamiltonian of the form 



_ 1 
oy ~ 2 



dx ^ (x) {—ivcr^dx + mcr2) {x) , 



where ai are the usual Pauli matrices and '5 {x) is a Majorana 2-spinor that writes 

( Ui^) 



(25) 



(26) 



with boundary condition ^^(0) = C-l(O)- The Hamiltonian ( p5|) is exactly-solvable being quadratic in terms of the 
fermions and the resulting eigenvectors read as follows in the Heisenberg representation 



^{x,t) 



1 



fe>0 



cos {kx + 9k) + i sin (kx) 
''^ ' cos{kx + 6k) — isin{kx) 



+ H.C. 



1 



'{m)r], 



(27) 



where is a fermion annihilation operator with k — irn/L, 77 is a zero mode Majorana fermion, and 6 is the Heaviside 
step function. In Eq. (27), denotes the energy dispersion of the model: 



Cfc = \/ u2fc2 _|_ ,^2^ 



(28) 



and 9k is given by 



cos 9k 
sin 9k 



m 



(29) 



For a positive mass m, one observes from the decomposition ( p7[ ) the existence of an exponentially localized state with 
zero energy inside the gap. Such localized Majorana fermionic states have already bie¥)3-,discussed in several different 
contexts such as^theJaolon edge state in an attractive one-dimensional electron gas,l 
fermion model,E3'C§EZI and the problem of a magnetic impurity in a superconductor 
out recently that such bound states may find applications in quantum computation 

^From this analysis of the toy model (|5|), we deduce the decomposition of the triplet and singlet Majorana fields in 
the basis of the eigenvectors of the Hamiltonian ( p3| ) subject to the boundary conditions (22). For Jj_ < i.e. nit > 
and TTis < 0, one obtains for the triplet sector a = 1,2,3 with obvious notations 




the random mass Majorana 
^' Finally it has been pointed 



(x,t) 



1 



ek 



cos {kx + 9l) + i sin (kx) 
cos {kx + 9\.) — i sin {kx) 



H.c. 




whereas the decomposition for the singlet excitations with m<; < reads 



{x,t) = 



1 



E 

fc>0 



COS {kx 
cos {kx 



9f,) + i sin {kx) 
9f.) ~ i sin {kx) 



+ H.C. 



(30) 



(31) 



Therefore, the localized Majorana zero mode state only appears in the triplet sector for a ferromagnetic interchain 
interaction. The translation of these results to the context of the spin-1 chain is straightforward: we just need to 
consider only the triplet sector. We find that in the Haldane gap phase {nit > 0) we have localized Majorana fermion 
modes at the edge, but, in contrast, the dimerized phase {rrit < 0) is characterized by the absence of such degrees of 
freedom. In Ref. |3^, it was shown that local zero modes were associated with kinks and antikinks of the mass m{x). 
In the present problem, with a semi-infinite system, the edge can be seen as a mass kink m9{x) and local zero modes 
should thus be induced irrespective of the sign of m. However, with the semi-infinite chain, the boundary condition 
( p^ ) selects only one chiral component. The sign of the mass then determines whether the local mode belongs to the 
physical chiral component. This is the reason for the difference of physical behavior between positive and negative 
masa^. 
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B. The uniform component of the magnetization profile 



With all these results, the physical properties of the edge states of the open spin-1 chain can be investigated. We 
first analyse the smooth part of the magnetization profile of the system. To this end, we consider the uniform part 
M(a;) of the total spin density (S+(x) = Si(a;) + 82(2;)) which takes the following form in the continuum limit using 
Eq. ®: 

2 

M(a;) =^(J„^(a;)+JaL(:E)). (32) 

a=l 

With the help of the identification (|l9|), we immediately find that the field M(a;) is expressed locally in terms of the 
Majorana fermions that account for the triplet excitations in the system: 

ix) = -'-e'^"^ en (x) en (x) - e'^'^^L (^) ^ (x) • (33) 
Using the decomposition (|30|), we write the uniform density M(a;) in the basis of the eigenvectors of the Hamiltonian 

nt (PI) 



iVr (x) = -ze^SyV— e"^"*"^/"* - W — e"™'^/"' V e'^''" cos (fca; + 6'f.) UArj' + H.c.) 
vt V ''^tL ^-^ 



k>0 



^abc 



- ^ E 1' ^) + ^ 5' ^) '^k^, + H.c) , (34) 

fc,g>0 



with 



A (fc, q, x) = cos [kx + cos (qx + 0*) + sin {kx) sin [qx) 

B (fc, g, x) — cos [kx + 0^) cos [qx + 6**) — sin {kx) sin (gx) . (35) 
The total uniform magnetization So is defined by 

/>oo 

So = / dxM {x) , (36) 

so that we obtain 







fc>0 

The first term in this-|equation describes a spin- 1/2 moment since it corresponds to the Majorana representation of 
a spin-1/2 operator.E3 In particular, the result ( |37| ) implies that the electron-spin-resonance (ESR) reponse of the 
cut two-leg spin ladder with a ferromagnetic interchain coupling decomposes into the bulk response and the response 
at a chain end. The latter one is identical to the ESR response of anjaojated spin-1/2 impurity. Since there is a 
continuity between the weak and strong coupling limits in this system,E3u3'Lj the Majorana approach provides thus 
an alternative description of the chain-eiid S=l/2 mode of the open spin-1 Heisenberg chain that has been obtained 
within the Schwinger bosons formalism. E3 

The uniform part of the magnetization profile of the model can be also read from the decomposition (^). For 
completeness, we give in the Appendix B an alternative derivation of the z-component of the uniform magnetization 
profile of the cut two-leg spin ladder without using the Majorana fermions method. We obtain the following result 
using Eq. (||) 

(M" (x)) = ^e-2"*^M \^ (38) 

Vt 2 

which can be interpreted as a spin-1/2 chain-boundary excitation localized over a length vt/2mt with an amplitude 
2mt/vt- This implies that the size of the spin-1/2 edge state diverges while its amplitude vanishes as the SU(2)2 
WZNW critical point of the S=l biquadratic chain is approached from the Haldane phase, in full agreement with 
the DMRG analysis of Ref. p5|. The Majorana fermion description also implies that the uniform component of the 
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magnetization profile should not be affected by temperature in the absence of an applied magnetic field. Let us finally 
mention that if the triplet mass mt is negative then the S=l/2 chain-boundary excitations disappear as it can be 
easily seen from the decomposition (|27|). We thus conclude that these free S=l/2 end spins are absent in a ladder with 
an antifcrromagnetic interchain exchange J_l > as well as in the dimerized phase of the spin-1 biquadratic chain. It 
is worth noting that the absence of fcee spin-1/2 moments in the spontaneously dimerized phase of a frustrated spin 
chain has been shown very recently.Ej 



C. Calculation of the NMR relaxation rate 



The NMR relaxation rate 1 /Ti of the cut two-leg spin ladder with a ferromagnetic interchain interaction can be 
computed by means of the Majorana approach described in the previous sections. For the standard two-leg spin 
ladder, it has been theoretically investigated in Refs. p5||5^ , |57| . We shall only consider here the uniform part of the 
NMR relaxation rate and restrict for simplicity to the contribution that identifies to the 1/Ti of the spin-1 Heisenberg 
chain in the limit of strong ferromagnetic interchain coupling —J±^Ji 



The general formula giving this NMR relaxation rate reads as follow: 

1 T 
—— = -Imx(x,w), 

ll[X} UJ 



(39) 



where uj is the nuclear resonance frequency which is the smallest energy scale of the problem: uj <C T,mt. We 
introduce the following susceptibility to perform the calculation of the NMR rate 1/Ti: 



iujn)^ / rfTe*"'"^ (T^M (x, r) • M (x, 0)) , 

"'0 



(40) 



with the analytical continuation x(a;,a;) — x(a;, zti;„)|ia)„->i^+jo- Using the decomposition (34) in the basis of the 
eigenvectors of the Hamiltonian Tit ( p3| ) that describes the triplet degrees of freedom, the NMR relaxation rate can 
be expressed as 



1 



QTtt 



q,x) [np (ffcj 



np {el)) Si 



(41) 



fe,9>0 



np{e) being the Fermi distribution function. The sum in Eq. (|4l] ) can be replaced by an integral through the 
substitution X]fc>o ^ ^ lo '^^/tt and the NMR relaxation rate simplifies, in the low-temperature limit T ^ mt, as 



Ti(x) 



A / dk e-^'^/^ 



T^Vt JO 



el 



[cos^ (kx + 9l) + sin^ {kx)f, 



(42) 



where the frequency uj insures the convergence of the integral at A: = 0. Using the energy dispersion ( p8| ) of massive 
Majorana fermions and the identification (E3) , one finally obtains the expression 



1 



Ti{x) 



TTV. 



t Jo 



ivtk)' 



2mtVtk 



+ 



mtVtk 



(vtk) + ml 



■ sin {Akx) 



{{vtkf 
2m1 



{vtkY 



(1 - cos(2fc.x)) + 



{vtk) 
ml {^{vtk) 



sin {2kx) 



2 {{vtkf +m 



(1 - cos {4kx)) 



(43) 



At the extremity of the chain (x = 0), the NMR relaxation rate takes a simple form in the low-temperature limit 



1 



6mt 



Ti {x = 0) TTwf 



T 
mt 



1 e 



T 
mt 



-mt/T 



(44) 



El (x) being the exponential integral function. Therefore, we conclude that the presence of the boundary leads to 
a narrowing of NMR linfi-.at low temperature compared to the bulk system. Ja principle, this NMR rate can be 
measured experimentally^ by measurements of nuclear magnetization recovery.Ej 
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Now, we turn to the calculation of the x dependence of 1/Ti. The sine terms that appear in Eq. (^3|) can be 
rewritten as 



12to, 



2mtX 



sinhfl 



6mt 



-•t JO 



In the regime T <ti mt, this expression can be approximated as 



where the function (p{y) is defined by 



cosh-" e 



■ sm 



sinhf 



2if 



(45) 



(46) 



n=0 



(271+1)!' 



,2n+l 



(47) 



and is a thermal length which writes 



6 



(48) 



This length scale diverges when T ^ and plays the role of an effective coherence length for the NMR relaxation 
rate. Similarily, the cosine terms of Eq. (km) can be rewritten in the low-temperature limit as 



de 



-ratS jlT 



3 — 4 cos 



6* + cos 



—( 



(49) 



We note that, for x ^ a i.e. far from the chain end, the low-temperature behavior of the NMR relaxation reads as 
follows 



1 



Ti (x > a) TTwf 



d0 



-mte^/2T 



.9^ + ^ 



(50) 



which corresponds to the bulk behavior of the NMR relaxation rate of the spin-1 Heisenberg chain found in Ref. ^ 
where the Haldane gap identifies to the triplet mass mt- 



IV. STAGGERED MAGNETIZATION AND DIMERIZATION PROFILES 



The staggered magnetization component of a two-leg spin ladder with a defect has been investigated semiclassically 
in Ref. |6l|. Such semiclassical approach has the inconvenient of breaking SU(2) rotational symmetry. Nevertheless, 
it gives useful qualitative indications on the expected magnetization profile. For the open ladder, the boundary 
condition on the bosonic fields is $±(0) 0. In the bulk, a semiclassical minimization of the ground state energy 
implies (<i>-|-) = {J± > 0) and ($+) = n/2 {J± < 0). Thus, we expect no staggered magnetization profile in the case 
of an antiferromagnetic rung coupling, and a profile with exponential decay far from the boundary in the case of a 
ferromagnetic rung coupling. In this section, we present an approach that has the advantage over the semiclassical 
method of preserving the full rotational symmetry. As is well known, the low-energy j^roperties of the two-leg spin 
ladder can be described using four decoupled off-critical two-dimensional Ising models. l3 In particular, this approach 
allows the calculation of the leading asymptotics of the staggered part of the spin-spin correlation functions which 
involve non-local operators in terms of the underlying Majorana fermions. In this section, we shall exploit the 
existence of a similar mapping for the semi-infinite two-leg spin ladder to determine the staggered component of the 
magnetization profile and the induced dimerization in the system. 



A. Staggered magnetization 



Let us discuss more precisely this mapping onto an effective Ising model. It is well known— hhat-a—lD theory of 
massive Majorana fermions describes the long-distance properties of ID quantum Ising model.E3Ho'EJ For a recent 
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detailed review on this correspondence, the reader may consult for instance the chapter 12 of the book Ref. |2^. In the 
case of a semi-infinite system, this inapping remains valid and the boundary conditions on the Ising spins depends on 
the ones for the Majorana fermions.E3 More specifically, we shall follow here the conventions of Ref. |6^ so that if the 
Majorana fermions £,r,l obey the boundary condition 

efl(o) = eL(o), (51) 

then the Ising model satisfies to a free boundary condition (i.e. the boundary spin is free to fluctuate and takes the 
values ±1). On the other hand, the Ising model experiences a fixed boundary condition (i.e. the boundary spin is 
fixed to the value (t(0) = 1 for instance) when the Majorana fields verify 

ei?(0) = -a(0). (52) 

The mass m of these fermions is a linear measurement of the deviation of the temperature with respect to the critical 
one: m = Tc — T as in Ref. p€ such that a positive mass corresponds to the low-temperature phase of the Ising model. 
The low-energy Hamiltonian ( 23|) of the cut two- leg spin ladder with the boundary conditions (|2^ ) on the fermions can 
thus be viewed as four decoupled off-critical ID quantum Ising models with free boundary conditions. In particular, 
the localized Majorana fermionic states with zero energy in the triplet sector found for a ferromagnetic interchain 
coupling ( J_L < 0) in Section III can be interpreted physically, in the Ising mapping, as a domain wall attached to the 
boundary which separates two domains of opposite magnetization (rrit = T^ — T > 0). In the singlet sector, one has 
in contrast rus < so that the corresponding Ising model with free boundary conditions is in its disordered phase. 
As a consequence, the zero-energy Majorana mode cannot exist in that case as it can be seen from the decomposition 

®- my 

The next step of the approach is to use the exact resultsE3EIle3 known for the semi-infinite Ising model to determine 
the staggered part of the magnetization profile of the cut two-leg spin ladder. To this end, the staggered magnetization 
n-|_ = rii + n2 of the total spin density S+ = Si + S2 is expressed in terms of the order and disorder operators o-frJ^ 
the different Ising models using the bosonic description (0) and the bosonization approach for two Ising modelsEjO 
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n_^_ ^ MiO'20-3Aio 

~ criCT2M3Aio- (53) 

At this point, it is worth discussing on the ground state degeneracy of the semi- infinite two-leg spin ladder with a 
ferromagnetic interchain coupling. As it was first pointed out by Kennedy,E3 an exponentially low-lying triplet, above 
the singlet ground state, is found in the Haldane gap for a finite open spin-1 Heisenberg chain. In the thermodynamic 
limit, the ground state is thus fourfold degenerate. At first sight, it seems difficult to reproduce this result starting 
from three decoupled semi-infinite quantum Ising models. Indeed, in the strong coupling limit —J±^ J||, the singlet 
degrees of freedom are frozen and the three Ising models for the triplet sector are all in their ordered phases (rrit > 
for J± < 0) so that (ui) 7^ (i = 1,2,3). In that case, each Ising model has a doubly degenerate ground state 
which gives thus an eightfold degeneracy. However, it is important to note that there is a redundancy in the Ising 
description since the triplet Hamiltonian in Eq. ( |23|), t he boundary conditions on the Majorana fermions (p2|), the 
Ising representation of the staggered magnetization (p3|) arc all invariant under the following transformation 



o-i -<Ji, (54) 

which leads to a physical fourfold ground state degeneracy as it should be. Let us return to the calculation of the 
magnetization profile for a ferromagnetic interchain coupling J± < where mt > and m^. < 0. The identification 
( ^ ) shows that the average staggered magnetization goes to zero far from the chain end since (/ii,2,3) = in the case 
of a positive triplet mass. However, due to the presence of the boundary, a staggered magnetization can appear close 
to the chain end i.e. when x = 0. The magnetization profile encodes the cross-over effect on the local magnetization 
as a function of the distance from the boundary. The magnetization profile of the spin-1 chain is obtained froni-4lie 
one of the ladder with ferromagnetic interchain interaction by taking the limit |r7is| ^ 00 or equivalently fiQ — > lE3'0 
In this respect, let us first present general results by exploiting the duality transformation on ID quantum Ising 
model. This transformation exchanges the order and disorder operators a ^ fj, but also the boundary conditions on 
Ising spins i.e. free boundary conditions become fixed and vice versa. Therefore, one obtains the following equivalences 
on the different one-point functions of the model 



{a (T > rc))froo = 



{fl (T < Te)) fixed = 
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(cr (T > Tc) (x))fixcd = {^J■{T < Tc) (a;))froc = cr^F 
(cr (T < Tc) (x))fixod = {tJ^{T > Tc) (a:))frcc = CooG 
(cr (T < Tc) (a;))frco = (M > Tc) (x))fixod = o-oc^i? (-^) 



(55) 



w being the velocity of the underlying Majorana fermion and CTqo is the expectation value of a (respectively /i) for 
T <Tc (respectively T_> Tc). An estimate of CToo valid for m < v/a is (Too = 2i/i2e~^/^^^/^(|TO|a/i;)^/^ where A is 
the Glaisher constant O It is indeed obvious that one has {(t{T > Tc)) hcc = for an Ising model with T > Tc and 
free boundary conditions. In contrast, one should observe that, even in the disordered phase of the model, a non-zero 
magnetization (cr(T > Tc))fixed 7^ exists for fixed boundary conditions since the Ising spins are polarized at the 
boundary. In that case, the precise cross-over between the boundary and bulk behaviors is described by the function 
F. The staggered part of the magnetization profile of the cut two-leg spin ladder with a ferromagnetic interchain 
coupling can thus be deduced from the correspondence (53) and the general results (|55|) 
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(56) 



which exhibits a full rotationally invariant form as it should be. Remarkably, a staggered magnetization appears 
although there is none for an isolated spin-1/2 Heisenberg chain. The expressions of (S-and H are exactly known and 
have been determined by BarievEll from a lattice description and later by Konik et al.E3 in the continuum case by the 
form factor approach. As it is shown in the Appendix C using this latter formalism, the function F can, in fact, be 
directly expressed in terms of G: 



F{x)^e-''G{x). 

As a consequence, the z-component of the staggered magnetization for instance simplifies as 



ml\ms\a*Y^^ _„^^/„^^2 f q f rrHx\ ^ ( \ms\x 



vfvs 



\ Vt 



Vt 



(57) 



(58) 



In the case of the open spin-1 chain, performing the substitution jiQ ^ 1 in Eq. (|53|), we obtain in a similar way the 
magnetization profile: 



{nl (x)) 



rnta\ ^-„rtx/vtjj2 f ^J}tx \ ^ f mtX 
Vt ) \ Vt ) \ Vt 



(59) 



form (see for instance 
Ea Complete expressions 



The functions G and H that appear in these equations can be cast into a Fredholm determi 
the Appendix C) or expressed in terms of a solution to the Painleve III differential equation, 
for G and H can be found in the Appendix For the sake of simplicity, we only need here the asymptotic behaviors 
of these functions which read as follows in the long-distance limit X = mx/v ^ 1E3E3 



G{X) 
H{X) 



-2X 



i60Fa:3/2 

1 ^-2X 



whereas in the short-distance limit X 



20FX1/2' 

mx/v <C 1, one has the following estimate 

G{X) - X-^/^ 
H{X) - X^/^. 



(60) 



(61) 
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^From these results, we deduce the behavior of the staggered component of the magnetization profile far from the 
chain end: 

(,;(^)) ^ (^!liKh!y^'e-«*^M ^ (^!^iKh!y^^-./6, (62) 

with a similar behavior for the open spin-1 chain. As expected, the local staggered magnetization decays exponentially 
with the distance from the boundary with a length scale that depends only on the bulk properties and identifies to 
the correlation length = vt/mt of the model. It is worth noting the absence of any x prefactor in front of the 
exponential term in Eq. ( |6^ ) in the long-distance limit a; ^ ^t- This suggests that the staggered magnetization for 
long chains with open boundary conditions is the relevant j-quantity to extract a very precise value of the Haldane 
gap as it has been done.^by means of the DMRG approach.Ea A similar exponential behavior is also obtained in the 
semiclassical treatmentEj and in a phenomenological theory of the open spin-1 chain describing the system as a spin- 
1/2 coupled to one-dimensional massive bosons.cj Comparing to Eq. (|3^), we notice that the staggered magnetization 
has a correlation length whereas lie uniform magnetization has a correlation length ^t/2, in agreement with the 
phenomenological free boson theory.E3 In contrast to the free boson theory, we find no x~^/'^ prefactor in the uniform 
component of the magnetization. In the short distance limit a; <C ^t, one obtains the following power law behavior 
from Eq. (^l|) for the two-leg ladder with a ferromagnetic interchain interaction: 

(n;(x))^^(aa;)i/2, (63) 
vt 

whereas, in the case of the spin-1 chain, this power law is modified to: 

Kx^^ - . 64 

We note that the staggered magnetization profile (^) has a vanishing intensity and a diverging correlation length 
when the Haldane gap goes to zero, in agreement with the DMRG analysis of Ref. Moreover, our calculation 
predicts that a staggerecLpagnetization will exist at T = in the — sector. This has indead been observed 
in a DMRG calculation. Ej At first sight, it seems to contradict the results of QMC simulations£3 However, these 
calculations are performed at finite temperature. For the one-dimensional quantum Ising model with free boundary 
conditions, there is no long range order in (cr) at T > due to the thermal nucleation of soliton excitations. Hence, 
we expect that as soon as the temperature is switched on, the average magnetization in S^^^ = state will vanish in 
agreement with what is observed in QMC calculations. 

The magnetization profile, in the antiferromagnetic interchain coupling case, can be investigated by a similar 
approach. For Jj^ > 0, one has now rrit < and iris > so that the Ising models of the triplet sector are in their 
disordered phases whereas the Ising model of the singlet sector belongs to its ordered phase. We thus obtain using the 
results (pSpS) that (n_|_) = and similarily it can also be shown that (n_) = (ni — = 0. Therefore, we conclude 
on the absence of S=l/2 chain-boundary excitations and of a non-zero magnetization profile for the cut two-leg spin 
ladder with an antiferromagnetic rung coupling. This result is consistent with the fact that the ground state of this 
model is always unique whether open or periodic boundary conditions are used. In this respect, the standard two-leg 
spin ladder with J± > 0, in contrast to the Jj_ < case, is not equivalent to the Haldane phase characterized by these 
S=l/2 chain-end excitatictfis eventhough they share similar properties like the presence of a spin gap, and a non-zero 
string order pacameter.EilO In fact, it has been recently pointed out that the two systems belong to two topologically 
distinct classes Jlj In particular, it has been argued that the S=l spin chain and the two-leg spin ladder with Jl > 
have two different types of string order that are intimately related to the valence bond structure of the ground states. 
The topological distinction is made by counting the number Qy of valence bond's crossing an arbitrary vertical line. 
In the case of an antiferromagnetic spin ladder, Qy is always even whereas it is odd for a system weakly connected to 
the spin-1 chain. Futhermore, the authors of Ref. ^ have noticed that for open boundary conditions ground states 
characterized by an odd value of Qy have spin- 1/2 edge states while these end states disappear when Qy is even. 
This is in full agreement with the results for the cut two-leg spin ladder obtained in this work within the bosonization 
approach. 
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B. Dimerization induced by open boundary condition 



The dimerization profile induced by the presence of a boundary can be also computed by this mapping onto semi- 
infinite Ising models. The dimerization operator in terms of the original lattice spins is defined by 

2 

e+n = (^1)" Sn,p • Sn+l,p- (65) 

p=l 

The bosonized description of this operator in the continuum limit reads as follows in terms of the bosonic fields ^± 
ofEq. (|) 

e+ COS cos (66) 



Using the bosonization representation of two Ising models J23C3'E3 this operator can then be expressed in terms of the 
different Ising disorder operators: 

e+ Afi//2Ai3Aio- (67) 

In the bulk, the system does not experience any dimerization pattern since the Ising models in the triplet sector are 
in their ordered phases for Jj^ < so that (/i;) = {i = 1,2,3). However, as for the existence-oLaJocal staggered 
magnetization, the presence of the boundary induces a non-trivial dimerization in the systemtZlEllLJ which can be 
obtained from the results (|55|): 

(,)) ^ ( !^iKk!) ,-3^..Mg3 (!I^)g( ^-^) . (68) 



Using the asymptotics (^, |6^) , we deduce the following estimates for the local dimerization 

V ^tVs J 

(e+(x)) ^ (9 (69) 

Note that the exponent 1/2 is identical to the exponent that would have been obtained in two decoupled gapless 
spin-1/2 chains by boundary conformal field theorytJ Physically, this means that the edge is making the system 
behaves as if it was gapless for distances shorter than the correlation length. In the case of the spin-1 chain, the 
dimerization operator ( |67[ ) simplifies to e+ ~ (/^o ^ 1) so that we get 

(e+(x)) ~ (^Y'e'^rn.-h^G' . (70) 



The long-distance limit of this dimerization has a similar form as in Eq. (|69|) and the short distance behavior is-modified 
to (e+ {x)) ^ x^^/^ . Again, this exponent could have been predicted from boundary conformal field theoryO We also 
observe that this exponent has been obtained in the DMRG study of biquadratic spin-1 chain at the SU(2)2 WZNW 
critical point .Efl 

In the antiferromagnetic interchain coupling case, a similar calculation can be made. The dimerization operator 
(H) has again a zero ground-state expectation value in the bulk since the Ising model in the singlet sector is in its 
ordered phase (m^ > 0). As seen above, the two- leg spin ladder with an antiferromagnetic rung coupling has no 
magnetic S=l/2 chain-boundary excitations but a localized Majorana state in the singlet sector still remains as it 
can be deduced from the decomposition ( p7[ ) with > 0. This zero mode manifests itself in the existence of a 
dimerization profile which is given by 



with the following asymptotics (^s — Vs/nis) 



e+ X - ^ 4 e =^/«%a;>C. 

V ""t^s ) 

(e+(x)) ^ {^-y\ (72) 
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In the case of the spontaneously dimerized spin-1 chain, the dimerization profile takes the form 

(,^(,))^f^y''G^f^V (73) 



Vt J \ Vt 



Its short distance asymptotics becomes thus: (e+ (x)) ^ x whereas the long distance one reads as follows using 
Eq. © 

(e+(x))^eoo 1+ ,;V h (74) 



16y^ x^/^ 



Coo being the non-zero bulk dimerization. 



V. EFFECT OF A STRONG EXTERNAL BOUNDARY MAGNETIC FIELD 

The effect of a strong applied magnetic field which fixes the spins at the boundary can be investigated using the 
Ising representation described in the previous section. To this endpJjCt us first recall the effect of a transverse edge 
magnetic field in the semi-infinite XXZ spin-1/2 Heisenbcrg chain.lZ3 It has been found that the system, along the 
entire XXZ critical line, renormalizes to the infinite field fixed point where the spin at the edge is polarized. In the 
bosonization language, one has an example of a c=l boundary fiow between the Dirichlet and Neumann boundary 
conditions. At the SU(2) invariant point, tli8-edge field is exactly marginal and a line of fixed point occurs between 
the Dirichlet and Neumann limiting casesJl3lI3 In the following, we shall only consider the physical situation where 
the spin at the edge is fully polarized or fixed so that it corresponds to the infinite field fixed point or Neumann 
boundary condition on the bosonic field <I>p associated to the spin-1/2 chain with index p = 1,2: 

a.$p(0,i) = OVi, (75) 

or equivalently it can be interpreted as a Dirichlet boundary condition on the dual field O^: 

Op (0, i) = Vt. (76) 

The value of the constant in this expression stems from the fact that a magnetic field along the x-axis is considered 
in the following. The actual direction of the applied field is not important since the model is SU(2) invariant. The 
chiral fields ^'±i?,L, defined by Eq. (H), are no longer independent due to the boundary condition (jT^) and satisfy now 

$±i (0) = $±,^ (0) , (77) 

from which we deduce the following analytic continuation {x > 0) 

$±i(x,t) = $±i^(-a;,^). (78) 



The change of boundary conditions in comparison to the Dirichlet case in zero field has several consequences. 
First of all, the commutator between the left and right bosonic fields is modified due to the folding condition ([T^): 
[^±r{x), $±l(2/)] — +iTT. As a consequence, the low-energy Hamiltonian of the model for a ferromagnetic interchain 
coupling Jj^ < is still given by Eq. ( ]2^ ) but with a negative triplet mass mt = J±X^/2t: < and a positive singlet 
mass rris = —'iJ\}?l2'K > 0. Moreover, the boundary conditions on the Majorana fermions have also changed in the 
Neumann case (^7|). They can be deduced as in Section II from the identifications (|l^; |l^) so that we obtain the 
following boundary conditions 

^i(O) = -eniO) 

?i(o) = ^(0) 

^i(O) = eniO) 

^LiO) = -aiO). (79) 

One can interpret these results in light of the Ising description presented in the previous section. The Ising models 
in the triplet sector with index 2, 3 (respectively 1) have free (respectively fixed) boundary conditions and belong 
to their disordered phases (mt < 0). The Ising model that accounts for the singlet excitations has fixed boundary 
conditions and is in its ordered phase (rn^ > 0). The S=l/2 chain-boundary excitations of the open spin-1 chain have 
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thus disappeared in the presence of a strong applied edge magnetic field. It remains only a single localized Majorana 
fermionic state with zero energy that describes fluctuations in the = triplet subspace which is unaffected by 



the applied magnetic field along the x-direction. The Ising representations of the staggered magnetization (53) and 



the dimcrization operator (67) have to be modified slightly due to the change of sign of the commutator between the 



left and right components of the bosonic fields $± and they are now given by 



£4 



(TiCr2CT3Cro- (80) 



^From these results and the identification ( p5[ ) , we thus deduce the staggered magnetization and dimcrization profiles 
of the two-leg spin ladder with a ferromagnetic interchain coupling in a strong applied magnetic field along the x-axis: 

{n% {x)) ~ e-l™'l-/^'G3 (^^^ G (^^) , (81) 

whereas (nV) = {n^) = = 0. The asymptotics of the non-zero staggered magnetization can be extracted from 
Eqs. 0) 

{n%ix)) ^ x~^'\ (82) 

The short distance exponent is identical to the one predicted-from boundary conformal field theory in a gapless 
spin- 1/2 chain with a strong magnetic field at the boundary.O The same result ( p^ ) also holds in the case of the 
spin-1 chain, albeit with a different short distance behavior {rv^ (x)) ^ x~^^^ which can be obtained from boundary 
conformal field theory. We conclude that the staggered magnetization in a strong applied field decays in the same way 
as in Eq. (|6^) far from the boundary but is enhanced in the vicinity of the chain end in comparison to the behavior 
( |63| ) in zero field. These repilts are in agreement with QMC simulations of the spin-1 Heisenberg chain with free and 
fixed boundary conditions 

A similar calculation can be made in the case of an antiferromagnetic interchain interaction. The only difference is 
that we must make the following substitution T — Tc ^ Tc — T. For a strong applied field along the x-direction, we 
get now 



ms\x 



(83) 

and (n^) = {n'\) = 0. However, there is now a non-zero staggered relative magnetization (n_ = rii — in the 
Jj_ > case which can be determined using the Ising representation of this operator: 



<Tia2fJ,3fJ'0, (84) 



so that 



and (n^) — 0. Finally, the dimcrization profile in the J± > case reads as follows 



(n^J^ (x)) ^ e-™'-/"'G2 ( H ( ^) H ( , (85) 



x^^/K^, (86) 



Vt J \ Vs 

SO that we obtain the following asymptotics respectively in the long and short distance limits 

(e+(x)) ~ e--/«= 



(e+(x)) ~ x'/\ (87) 



16 



We close this section by discussing the case of the spontaneously dimerized spin-1 chain. The Ising representations 
of staggered and dimerization fields are now given by Eq. (|80| ) with (Tq — > 1. The staggered magnetization and 
dimerization profiles are 



In the short distance limit, we obtain the following power law behaviors 



whereas for long distances we have 



{e+{x)) ^ x'/\ (89) 
(nlix)) - e-2^/«' 

/ .1/2 2x/6\ 

^ e^{l-^-^^T7^j. (90) 

Therefore, we observe that the dimerization reaches here its bulk expectation value from below in contrast to the free 
boundary case (u4)- 



VI. CONCLUDING REMARKS 



In this paper, we have investigated the nature of the chain-boundary excitations of the cut two-leg spin ladder|-a|ad| 
the open spin-1 chain by means of the bosonization method. The crucial point of the analysis is the mapping23'0 
of the low-energy Hamiltonian of these systems onto free massive Majorana fermions (or equivalen t ly , d ecoupled 
non-critical quantum Ising models) with suitable boundary conditions. In particular, the exact result^r^r tE3 of the 
semi-infinite one-dimensional quantum Ising model allow the determination of the low-energy properties of the cut 
two-leg spin ladder such as, for instance, the magnetization and dimerization profiles. For a ferromagnetic rung 
coupling {J± < 0), the system is characterized by the presence of fractional spin-1/2 edge states which, in the limit 
J±/J\\ ^ — oo, identify to the well known S=l/2 chain-end degrees of freedom of the open spin-1 chain. In this 
respect, the approach, presented in this paper, pro\pies an alternative derivation of the existence of |these edge states 
first predicted theoretically within the VBS modelt3 and the Schwinger boson mean-field analysisO In the case of 
an antiferromagnetic interchain interaction Jl > 0, no S=l/2 chain-end excitations are found but a non-magnetic 
localized Majorana fermion zero mode is still present and leads to the formation of a non-zero dimerization profile in 
the system. 

The magnetization and dimerization profiles, derived in this paper, should be confronted to numerical simulations 
of the cut two-leg spin ladder or the open biquadratic spin-1 chain in the vicinity of the SU(2)2 WZNW critical point. 
Due to the semi-infinite geometry considered here, our results would best be compared with DMRG .calculations of a 
finite spin-1 chain with a spin-1/2 attached to one of the extremities to cancel one of the edge states.llBc3 At this point, 
it is worth noting that all the calculations were done at zero temperature. Our results could in principle be extended, 
to finite temperature using the thermal form factor techniques derived for the one-dimensional quantum Ising model.O 
Unfortunately, it is not an easy task within this formalism to obtain explicit expressions for {(j{x)) for fixed boundary 
conditions. This makes difficult any direct comparison to QMC simulations.EScJ However, one can argue by finite size 
scaling arguments that the correlation functions should not be strongly affected by a finite temperature as long as 
TO ^ T. This can be checked by an explicit computation of the two-point correlation function.ll3 We also stress that 
at finite temperature and for free boundary conditions, one has {(t{x)) = in the quantum Ising model. This implies 
the absence of any staggered magnetization profile in the S — A state and also of a non-zero string order parameter 
at finite temperature in agreement with the QMC simulations.li3 

Regarding perspectives, the approach, presented in this work, could be applied to other one-dimensional gapful 
systems. The effects of an uniaxial single-ion anisotropy ^^iSf)"^ on the magnetic properties of the open spin-1 
chain can be investigated. Since the different species of Majorana fermions do not-interact, we expect that spin-1/2 
edge states excitations should still be observed, in agreement with the QMC results.LJ The calculation of magnetization 
profiles with our method should not pose any difficulty. j-The approach could also be generalized to study the effect 
of a weak bond or magnetic impurities in a spin-1 chainEII as it will be discussed in a separate publication. Another 
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interesting situation is-the nature of the edge states of two open spin-1/2 Heisenberg chains coupled by a biquadratic 
interchain interaction.ES Due to the extended 0(4) symmetry of the model, we expect to find two spin-1/2 excitations 
at the edge of the system. A more challenging problem is the generalization of our approach to S* > 3/2 Heisenberg 
chains. According to Ref. |l^, it is expected that edge states with fractionalized spin S" exist in the spin-S chain 
with S' — S/2 (respectively S' = {S — l/2)/2) for integer (respectively half-odd-integer) spins. This conjecture, 
based on the large-N limit of SU(f^ quantum antiferromagnets and strong-coupling expansion, has been verified by 
a DMRG analysis for S = 3/2, 2E3 A generalization of the approach presented_^here in the S=l case is to describe 
spin-S Heisenberg chain as perturbed SU{2)2s = U{1) ® Z2S WZNW models.EZlE^ For S half-odd integer, only the 
parafermion sector Z2S is gapped. We should thus expect the edge spin excitations to be generated by the bound 
states of the massive Z2S theory on the half-line. For S integer, edge excitations should be induced by the boundary 
bound states of the perturbed WZNW model. A similar problematic should also be considered for the related problem 
of n-leg spin-1/2 ladders. Finally, an interesting question would be the study of interactions between edge states in 
chains of finite lengthE3 using the Majorana fermions description. 
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APPENDIX A: BOSONIZATION APPROACH OF THE OPEN S=l/2 HEISENBERG CHAIN 

In this ap j pej Bdia, we describe the bosonization approach of the spin-1/2 Heisenberg chain with open boundary 
conditions .Eacj'Ea'c3 This enables us to fix the conventions that will be used in this paper and also to discuss some 
subtleties related to the presence of open boundary conditions. To this end, we consider the repulsive Hubbard model 
at half-filling with open boundary conditions described by the Hamiltonian 

N-l N 

Hu^-tY, (cLc»+i<T + i?.c.) + C/ ^ Tii^n.i , (Al) 



where Cia is the electronic annihilation operator of spin index a at site z (1 < i < N) and rii^ — cl^Cia 

stands for the occupation number of electron with spin index a. The summation over repeated greek symbols is 
assumed in the following and the hopping term t is positive. In this model, it is well known that a charge gap rric 
exists for any positive value of the interaction U and in the low-energy limit {E <C rUc) only the spin excitations 
remain and describe the universal scaling properties of the spin-1/2 Heisenberg chain. In this way, we shall derive the 
continous description of t he spin density of the S=l/2 Heisenberg chain with open boundary conditions starting from 



the electronic model (Al). An alternative approach as described in Refs. |35| , p8| , [l0| is to consider the spin-1/2 XXZ 
Heisenberg chain with open boundary conditions and the use of the Jordan- Wigner transformation. Since in this work 
we shall only consider SU(2) invariant interactions, it is more appropriate to start from the Hubbard mo del ([Al]). 
The open boundary conditions are taken into account by introducing two fictious sites-|Q-and A'^ -|- 1 in Eq. ( |Al[ ) and 
by imposing vanishing boundary conditions on the fermion operators: cq = pf-i-i = OEj'EZI The low-energy properties 
of the model can ±lieB|-he--determined by applying the bosonization methods with suitable boundary conditions on 
the bosonic fields .E§E30cj 

1. Non-interacting case 



In the low-energy limit, the continuum version of the non-interacting part of the Hamiltonian (Al) can be derived 



by expressing the lattice fermions Cna in terms of left and right moving spinful fermionic fields '^L,Raix)'. 

^ ~ i^^^ (x) + {-if/^ -^L. {x) , (A2) 

with X ~ na, a being the lattice spacing. The resulting boundary conditions on the fermionic fields of Eq. ( [A2| ) are 
thus 

^lAL) = -i-lt^^'^RAL), (A3) 
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with L = (TV + l)a. The left and right excitations are no longer independent due to the presence of these boundaries. 
The next step of the approach is the introduction of right and left moving bosonic fields ^r.lct through 



ITTTcr /4 



'27ra 

pi7rT„/4 



La 



/27ra 



(A4) 



where Ko- are Klein factors that obey the anticommutation relations {k^, ^o-/} = ^a-.ar to ensure the anticommutation 
between the fermion fields of different spin index. In Eq. (A4), we have also introduced some phase factors with 
T| = 1 an d T i = — 1 for later convenience. The boundary conditions on the chiral bosonic fields are then obtained 
fromEq. Q: 



(0) 



^^R, (0) + TT 



L 



(L) + TT - - 1 + 2q„n, 



(A5) 



Qa- being an integer. In our conventions, the total bosonic field ^„ with spin index cr and its dual are related to 
the chiral components through 



1 



Rcr 



(A6) 



so that Eq. ( |A5| ) imposes Dirichlet boundary conditions on the bosonic field: 

(0) 



TT 

2 



(A7) 



The low-energy dynamics of the non- interacting Hamiltonian TLq of the origin al m odel (Al) is thus described by 
two independent free massless boson Hamiltonian with the boundary conditions (A7): 



dx 



(AS) 



vp being the Fermi velocity and W„ is the momentum operator conjugate to ^o-- The mode decomposition of the 
bosonic field ^„ compatible with these boundary conditions reads as follows 



TT f L 



'TT TTOc 



X 

L 



(A9) 



where fc„ — mr/L, ana is the boson annihilation operator obeying [ana, ctlna/] — Sn,mSa.ai and the zero mode operator 
ttq^ has a discrete spectrum ^Jixq^. The mode deco mp osition of the momentum operator Ho. = dt^ajixvp conjugate 
to the bosonic field can thus be deduced from Eq. (A9) 



Her (x, = X! "T~ ^^"^ 



(AlO) 



In particular, one can check that the mode decompositions (A9, AlO) satisfy the canonical commutation relation: 

(x, t) , n^, (y, t)\ = i(5^,^,(5L(x - y), (All) 

— y) being the delta function at finite size: Sl{x) = J2n e''^"^/2i. The dual field Qa satisfies dxQa = ttIIo- and 
dfOa — vpdx^a so that one obtains the following mode expansion: 



E 



cos {knX) 



(A12) 
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where the zero mode coordinate 4'oa- is conjugate to ttq^: [^oo- , tToo-' ] ~ ^^aa' ^^'^ ^'^^ fixed by the boundary 
conditions ( [A7|). Finally, the mode decompositions of the chiral bosonic fields ^u^La can be determined by the 
identification ( [Aq ) : 



TT 






2 + 






TT 




(L 


2 + 


(1 


\ a 



n— 1 ^ 

oo 



L 

X — 



~ikn{x+VFt) _ t ifc„(a:+i>Ft) 



'T^ <P0o 



ikn{x-VFt) _ t -ik„{x-VFt) 



(A13) 



In addition, one can show that these chiral fields satisfy the following commutation relations when L 3> a 

[<^La[x,t) ,<^Lai[v,t)] = -iTTfJo-^cr/ sgu (a; - y) 

[(^R„ (x, t) , $_rct, (y, i)] = i7r(5cr,o-/ sgn (x - y) 

[^R„{x,t) ,^Lai{v.t)] = ifx = y = 

= —217: 5 a Ml ii X = y = L 

= -ZTTf^o^^c,, ifO<a;,2/<L, (A14) 

sgn(x) being the sign function. 



2. Effective spin density 

The next step of the approach is the introduction of the bosonic fields that describe the charge and spin degrees of 
freedom: 



LI 



V2 



V2 



(A15) 



This basis as well as the commutation relations ( A14 ) allow us to express the Hamiltonian ( [AS| ) in terms of two 
commuting gapless spin and charge contributions: 



Hq = 



VF 

27r 



C dx ({dx^cf + (ttH,)") +^£dx ((9,$.)' + (ttH,)') 



(A16) 



As well known, a weak Hubbard interaction preserves this famous spin-charge separation and opens at half-filling a 
mass gap rric for the charge degrees of freedom. In the spin sector, the effect of the interaction is exhausted by a 
renormalization of the spin velocity and by the existence of a marginal irrelevant contribution in the Hamiltonian. 
In particular, the interaction does not renormalize the bosonic field $s since it is protected by the underlying SU(2) 
symmetry of the model. Neglecting the logarithmic corrections introduced by the marginal irrelevant term, the 
low-energy (E <C mc) Hamiltonian that describes the universal properties of the S=l/2 Heisenberg chain is simply 



ns = ^ [ dx ({dx^sf + (vrH, 



(A17) 



Vs being the velocity of the spin collective mode. The boundary conditions of the bosonic field can be obtained 
from Eqs. (p^, |A15D: 



*.(0) = 



qn 

71' 



(A18) 



q being an integer. Similarily, the mode decompositions of the chiral bosonic fields ^sR,l read as follows with help of 
Eq. ([aT3|) 



oo . 

n— 1 ^ 

oo 



L 

X — Vst 



ikn(x—Vst) I —ikn{x — Vst) 



(A19) 
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with [(j)Qs,TTOs] 

relations 



i and [a„s , ct. 



t 1 

msJ 



: Sn^m- Moreover, we deduce from Eq. ( A14) that these fields obey the commutation 

[^sL{x,t) ,<i>sL{y,t)] = -insgn{x-y) 

[<i>,R{x,t) ,<i>,H{y,t)] = iT:sgn{x-y) 

['^sR{x,t),<^,L{y,t)] = ifa; = y = 

= —2iTT a X ^ y = L 

= -in if < x,y < L. (A20) 

At this point, it is important to note that the last commutator in Eq. ( |A2(]| ) has the opposite sign of the prescription 
made in Refs. |2^:^- The actual value of its sign stems from the fact that we are considering open boundary conditions 
in the lattice system which identify to the Dirichlet boundary conditions (A18) on the bosonic field ^g- In fact, one 
can derive the value of the commutator [^sR, ^sl] by a different method. The boundary condition at a: = on the 
chiral bosonic spin fields is 

$3L(0,i) = -$sfl(0,i) Vi. (A21) 
Since {x, t) = $sl {x + Vst) and {x, t) — ^sR {vst — x), one thus obtains the folding condition {x > 0): 

= (A22) 



which is satisfied by the mode decompositions ( |A19| ). Moreover, the commutator [^sR{x,t), ^sniy, t )] is fixed by the 
requirement that $s and Xlg are canonical conjugate operators so that by using the folding condition ( A22) we deduce: 



{x,t),^sL iy,t)] 



- (x, i) , i~y,t)] 
-in sgn (x + y) — —in. 



(A23) 



ft turns out that the sign of this commutator is important for the investigation of the S=l/2 chain-boundary excitations 
of the open two-leg spin ladder as described in Sections II and III of this work. 

With all these results at hands, it is straightforward to derive the continuum description of the spin density starting 
from the lattice spin operator S^: 



S: = I cKpC^p. (A24) 

a'^ia — X, y, z) being the Pauli matrices. Using the decomposition (|A2| ), the spin density separates into a uniform and 
staggered parts in the continuum limit: 



S {x) =3sR (x) + JsL (x) 



c/a 



Us (x) 



(A25) 



with the identification 



ja 



sR.L 



(A26) 



The bosonized description of the spin density can then be obtained with help of the bosonization formula ( A4), the 
commutation relations (A14), and the canonical transformation (A15). The resulting expressions for the uniform part 
read as follows 



•JsR 



1 



sL 



sR 



J, 



J, 



sR 

t 

sL 



2nV2 
1 

2nV2 



2na 
2na 



(A27) 
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whereas the staggered part is given by 









na 


nl = 


* r 




na 




A 


< = 


sm 




TTfl 



sm 



)s (V2 e, 
(V2 e,) 



(A28) 

A being a constant stemming from the charge degrees of freedom that have been integrated out in the low-energy 
regime E <C m^. The product HfHi has no dynamic and in this work we use the prescription — i for simpUcity. 



FinaUy, as it can b e checked from the correspondence (A27), the left-moving contribution of the uniform part of the 



spin density (A25) obeys the following operator product expansion (with a similar result for the right-moving term) 

i^abcjb 



(z) [w) 



' sL 



(w) 



8TT^{z-wf 2t: (z-w) 



(A29) 



with z ~ VsT + ix. The uniform left spin density Jsl identifies to the SU(2)i Kac-Moody currents which arc the 
generatofSjOf the conformal field theory associated to the criticality of the spin-1/2 Heisenberg chain (see for instance 
the bookca for a review). 



APPENDIX B: ALTERNATIVE DERIVATION OF THE UNIFORM MAGNETIZATION PROFILE 



In this Appendix, we derive the z-component of the uniform magnetization profile of the cut two-leg spin ladder 
with a ferromagnetic rung coupling without using the Majorana fermion formalism. To this end, we return to the 
complex fermion Hamiltonian ?i+ (17) with Jj^ < 0. The z-part of the total magnetization density is given by: 
=: ip\.j^ip+R + '4>\]^ip+L '■■ Using the boundary condition ( p^ ) and the Hamiltonian (17) with m = —J^j^-n > 0, 
we obtain the following mode decomposition for the fermionic fields ip+n^L- 



^+r{x) 



V 



= '^if{k,x)ak.+ + f*{k,x)ak-) 
k 

= ^{f*{k,x)ak^+ + f{k,x)ak^^), 



(Bl) 
(B2) 



where f{k,x) = cos{kx + Ok) + isin(fcx). Ok being defined by Eq. (p9[). The Hamiltonian 7Y+ ( |17| ) can be expressed 
in terms of the ak,± fermionic modes 



= ^ e(fc)(4_+afc,+ - 4__afc,_), 



(B3) 



with the energy dispersion e(fc) = \J v'^k'^ + m^. The uniform magnetization profile along the z-axis is then given by 



(M-(x)) = \y\co^\kx + Ok) + sin2(fcx) - l]((4.-afc.- + 4,+afc,+)) + — e-2™^/^(4ao). 



(B4) 



fe>0 



Using the expressions (|29|), noting that np{e{k)) + np{—e{k)) = 1 and J2k>o ^ ^ Io° dk/ir in the large L limit, the 
uniform magnetization simplifies as 



{M^{x)) = ^e-2-/^(4ao) 



V 

°° dk 





{vk)^ 



■ cos(2fca;) 



(vk)^ 



■ sin(2fca:) 



Performing the integrals, we finally find the following result 



2™ 

M^(a;) = e 

V 



-2mx/vl/ A 



[{alao) - 1/2]. 



(B5) 



(B6) 
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This profile is identical to the one obtained by the Majorana fermions calculation (p8|). 

If we apply a uniform magnetic field along the z-axis, the Hamiltonian 7i_|_ (B3) is modified as follows 

7^+ = ^ (re(fc) - h) : 4^^afc,, : -/i(aJao - 1/2), (B7) 

fe,r=± 

whereas the Hamiltonian ?i_ in Eq. (^) is not affected by the magnetic field. The resulting free energy per unit 
length is then given by 

= — V ln(l + e-'3(^('=)-'-'^)) - - ln(2 cosh{rjh/2)) = /^'^ + ff^^. (B8) 

Jo ^ P 

^From we recover the usual susceptibility and specific heat of the spin-1 chain. We see that f°'^fi° is the free 

energy of an isolated spin-1/2. Thus, for /i <§C /3 and h -^m, the thermodynamic properties of the systepi.are identical 
to those of an isolated spin-1/2. This result is in agreement with, the QMC simulations of long chainal3 and DMRG 
calculations of effective interaction of edge states in long chainsE3. 



APPENDIX C: CALCULATION OF THE FUNCTION F{x) 

In this Appendix, we compute the one-point function of the disorder operator in the low-temperature phase of 
the semi-infinite one-dimensional ouantum Ising model with free boundary conditions. To this end, the form factor 
approach to correlation functionsEJ will be used as in Ref. |6^ for the calculation of the magnetization one-point 
function. |— . ■_. 

Let us first recall some results on the form factors of the bulk quantum Ising modehO'Ej The excited states of this 
model are created by acting on the ground state with fermion creation operators A^: 

|0i...0„) - At((?i)...^t(0,^)|O), (CI) 

where the ^^'s are the usual rapidity variables parametrizing momentum and energy as p{Oi) = msinh^i, e{9i) — 
mcosh^i, m being the fermion's mass (m > for T < Tc) and its velocity has been set to unity here for simplicity. 
The creation and annihilation A operators satisfy the fermionic anticommutation relation normalized as follows 

{A((?i),At(02)} = 27rJ(0i-02). (C2) 
For T < Tc, the form factors of the order operator a are: 

(O|a(O)|ei...02„) n *^"^(^^)' (^3) 

l<i<j<2n ^ ^ 

whereas the form factors with an odd number of rapidities are zero. They are normalized such that the conformal 
limit of the spin-spin correlation function is 

(^W^(0)) = |^, (C4) 

with r = \/x'^ + t'^ and T = 'l~^l^'^e^l^ A~^l'^ra~^l^ , A being the Glaisher constant. In the low-temperature phase, 
the form factors of the disorder operator are given by 



(O|At(O)|0i...02«+i) n *^''^(^T^)' 

l<i<j<2n+l ^ ' 

and those with an even number of rapidities are zero. For T > Tc, the roles of a and [i are interchanged. 

With these results, one can extend the method of Rcf. ^ to calculate the one-point function of the disorder 
operator in the low-temperature phase of the semi-infinite Ising model with free boundary conditions. The free 
boundary condition on the Majorana fermions |(||l(0) = Cii(O)) is interpreted as a boundary state \B) which encodes 
all informations about the boundary condition£a In this approach, the Hilbert space of the theory is the same as in 
the bulk so that the one-point function can be extracted through 

(M(x))=5](0|A.(0)|n)(n|i?)e-^^", (C6) 
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\n > being a complete set of states of the bulk Hilbert space. The boundary state corresponding to the Ising model 
at T <Tc on the half line with free boundary conditions ia£3 



\B) = (1 + At (0)) exp 



|o>, 



(C7) 



with R{0) = — i coth((?/2). This boundary state contains a zero-momentum one-particle state which corresponds to a 
domain wall, attached to the boundary, that separates two domains of opposite magnetization (T < Tc). Such term 
contributes to the expectation value (C6) whi le it does not enter the calculation of the one-point function of the order 
operator. Expanding the exponential in Eq. (|C7|), we get the following expression using the fact that the form factors 
of /i are non-zero only for an odd number of rapidities: 

°° 1 roc in roc in 

(m(^)>=EA / ^(O|MO)|O;-0i,0i;...;-0„,0„)m)...i?(0„)e--[i+^™^'^('^^)^ 



from which we deduce the identity 



(C8) 



(C9) 



The next step of the approach is to use the form factor of fi (|C5| ) to derive an expression for A^mx). First of all, one 
has 



(0|At(0)|0; 01, -Oi; 9,,, -6l„) = J^tanh^ ^ tanh 6*, J]^ tanh^ 



2 



tanh'^ 



so that, using the expression of R{0), we obtain: 



cosh(Sfc) 



with 



2 cosh Qi cosh Qj 
cosh6'i -I- cosh 61, ■ 



Following Rcf. pSL we introduce the quantity 



V[p,,0,,mx)- cosh 0,-1- cosh 0j 
The function A(mx^ can then be expressed as a Fredholm determinant: 



(CIO) 



(Cll) 



(C12) 



(C13) 



°° 1 roo in roc in 

Aimx) = y A / — ••• / —detV{9,,e,,mx) 



V 

Det ( 1 + — 
Zir 



(C14) 



Using the results obtained in Ref. |68|, A{mx), given by Eq. (C14), coincides with the Fredholm determinant repre- 
sentation of the one-point function of the Ising magnetization for T < Tc with fixed boundary conditions i.e. G{mx) 
in our notations (see in particular Eq. (^5|)). Therefore, from Eq. (C9), we finally deduce the following relation: 

F{mx) = e^'"^G(ma;). (C15) 



APPENDIX D: EXPRESSION OF THE FUNCTIONS G{x) AND H{x) IN TERMS OF SOLUTIONS OF 

THE PAINLEVE HI EQUATION 

According to Bariev,0 the functions G and H, describing the cross-over effect on the local magnetization of the 
semi-infinite Ising model at T < Tc with free and fixed boundary conditions, can be expressed in terms of a solution 
r]{d) of the Painleve III differential equation. This latter equation reads as follows 

• 1] - —. (Dl) 



rjde^ \r]dej erjde 
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The boundary conditions on rj are: 



V{0) 



In -^+1E 



1 — [e oo), 

Zn 



(D2) 



jE being the Euler's constant. The functions G and H that are the building blocs of the staggered magnetization 
and dimerization profiles are related to the solution r]{6) by 



H{y) = ii^'\y)e^Y> 



1 



d0 <! ^v-'io) \ (1 - v'ieyr - ( ^ ) ) - -iv-'{e) - 1) 



dr/ 



1 



(D3) 



There is an interesting connection between the Painleve III differential equation and the two-dimensional s inh - 
Gordon equation. Indeed, the relation is obtained by considering ri{9) — e^^^^-* so that the differential equation (Dl) 
takes the form 



(fx Idx 
dO"^ 6 d9 



— 2 sinh2x. 



(D4) 



The functions G and H in Eq. (D3) can then be expressed in terms of x 



G{y) 



ex(»)/4gxp 



4 sinh X ^ 



4 sinh X 



dx 
dO 



(e>^W - 1) 



(D5) 
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